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^^ ■ Anyon systems are studied in connection with several interesting applications in- 

fT— , , eluding high Tq superconductivity and topological quantum computing. In this work 
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c/3 ' we show that these systems can be realized starting from directed polymers braided 
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together to form a nontrivial link configuration belonging to the topological class of 
plats. The statistical sum of a such plat is related here to the partition function of 
a two-component anyon gas. The constraints that preserve the topological configu- 
ration of the plat are imposed on the polymer trajectories using the so-called Gauss 



Cn . linking number, a topological invariant that has already been well studied in polymer 



physics. Due to these constraints, short-range forces act on the monomers or, equiv- 
alently, on the anyon quasiparticles in a way that closely resembles the appearance 



^T . of reaction forces in the constrained systems of classical mechanics. If the polymers 

^^ are homogeneous, the anyon system reaches a self-dual point, in which these forces 

vanish exactly. A class of self-dual solutions that minimize the energy of the anyons 
is derived. The two anyon gas discussed here obeys an abelian statistics, while for 
quantum computing it is known that nonabelian anyons are necessary. However, this 
is a limitation due to the use of the Gauss linking invariant to impose the topological 
constraints, which is a poor topological invariant and is thus unable to capture the 
nonabelian characteristics of the braided polymer chains. A more refined treatment 
of the topological constraints would require more sophisticated topological invariants, 
but so far their application to the statistical mechanics of linked polymers is an open 
problem. 



I. INTRODUCTION 



Knots and links are a fascinating subject and are researched in connection with many 



concrete apphcations both in physics and biology |l|-|25|. In this paper we study the statis- 
tical mechanics of a system of two entangled polymer rings. Mathematically, two or more 
entangled polymers form what is called a link. Single polymer rings form instead knots. 
We will restrict ourselves to systems in the configurations of 2s— plats. Roughly speaking, 
2s— plats are knots or links obtained by braiding together a set of 2s strings and connecting 
their ends pairwise [25^]. A physical realization of 2s— plats could be that of two rings topo- 
logically entangled together and with some of their points attached to two membranes or 
surfaces located at different heights. In nature 2s— plats occur for example in the DNA of 
living organisms jlll . 



27| . Indeed, it is believed that most knots and links formed by 
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DNA are in the class of 4— plats jjjj. These biological applications have inspired the research 
of Ref. |28| . in which 4— plats have been studied with the methods of statistical mechanics 



and field theory. In particular, in 



28| it has been established an analogy between poly- 



meric 4— plats and anyons, showing in this way the tight relations between two component 
systems of quasiparticles and the theory of knots and links. After the publication of 28] . 



interesting applications of analogous anyon systems to topological quantum computing has 
been proposed 129 



-l3l|. These applications are corroborated by the results of experiments 
concerning the detection of anyons obeying a nonabelian statistics, see for example |32| . 
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33|,_other systems 



While these results have appeared in 2005 and are still under debate 
in which non-abelian anyon statistics could be present have been discussed |34l . |35|. In the 
present case, the topology of the original two-polymer link is distinguished by the Gauss 
linking invariant, which can be obtrained from the amplitudes of an abelian BF model J36|. 
This implies that the statistics of the quasiparticles treated here is purely abelian. However, 
also abelian anyons may be exploited for quantum computations as it has been argued in 
Ref. Is 



Motivated by these recent advances, we study here the general case of 2s— plats formed 
by two polymer rings. Among all knot and link configurations, the class of 2s— plats is very 
special. For instance, it is possible to decompose the trajectory of a 2s— plat into a set of 
2s open subtrajectories that can be further interpreted as the trajectories of 2s polymer 
chains directed along a special direction. Without losing generality, we may suppose that 



this direction coincides with the z— axis. When the system of directed polymers is mapped 
into a field theory, a model describing anyon quasiparticles is obtained. The z coordinate can 
be related to "time", while the monomer densities of the 2s directed polymers become the 
quasiparticle densities of a multicomponent anyon model. A remarkable feature of polymers 
in 2s— plat configurations is that they admit self-dual solutions in which the energy of the 
system is minimized 



28] . Here we show that these solutions can be explicitly constructed 
by solving a sinh-Gordon equation. The conformations corresponding to such solutions 
should be particularly stable and thus observable, at least in principle. With the present 
technologies (38|, in fact, it is possible to realize polymer 2s— plats in the laboratory. 

Another advantage of restricting ourselves to 2s— plats configurations is that it is possible 
to distinguish their topological states in a more efficient way than what one could achieve in 
the general case of two linked polymer rings. Let us recall at this point that the trajectories of 
real polymers are impenetrable and thus, if no rupture occurs, they are bound to stay in the 
initial topological state while subjected to thermal fiuctuations. However, in the Edwards' 
model used here polymers are "phantom" [39'j. Without any control, their trajectories are 
allowed to cross themselves and thus the global topological configuration of the system may 
change. To find a powerful and reliable method in order to forbid such changes of topology 
is the most difficult problem of the statistical mechanics of polymer knots and links. Up to 
now there is no analytical model that is able to deal with the statistical mechanics of polymer 
knots. For this reason, in this work we assume that each polymer ring composing the link 
can be in any knot configuration. Only the topological configurations of the link formed by 
polymers belonging to the 2s— plat will be distinguished. This goal is achieved by using the 
Gauss linking number in order to impose the necessary topological constraints. The Gauss 
linking number is a topological invariant given in the form of a double contour integral, where 
the contours coincide with the polymer trajectories. Unfortunately, it is a weak topological 
invariant, so that many nonequivalent topological configurations characterized by the same 
value of the Gauss linking number are allowed. However, once we restrict ourselves to a 
given 2s— plat, we are implicitly imposing a much more stringent topological condition on 
the system. Indeed, its topological states are in this way not only limited by the value of 
the Gauss linking number, but are also forced to vary within the much smaller set of states 
that are compatible with the structure of the 2s— plat. 

Even if the Gauss linking number is one of the simplest topological invariants, its expres- 



sion is very complicated. As a consequence, after imposing the topological constraints, the 
action of a system of topologically entangled polymers becomes both nonlocal and nonpoly- 
nomial. The nonlocality is due to the double contour integral over the polymer trajectories. 
The nonpolynomiality arises from the fact that the integrand is a nonpolynomial function 
of the components of the radius vectors determining the positions of the monomers in the 
space. The situation is somewhat reminiscent to that of holomic constraints in the classical 
mechanics of particles. When these constraints are fixed by means of Lagrange multipliers, 
within the particle action new terms appear which are related to the reaction forces. The 
striking difference is that topological constraints are not holonomic and have "memory" 
in order to keep track of the global conformation of the chain. This last property causes 
the nonlocality of the action after fixing the constraints. The price to be paid to recover 
locality and to have a standard action is to introduce topological fields which interact with 
the monomers in such a way that the topological configuration of the link is preserved. To 
some extent, these interactions may be considered as the equivalents of the reaction forces in 
classical mechanics. The passage to the topological field theory description is not straight- 
forward. In particular, it requires to find a topological field theory with an amplitude of 
metric independent and gauge invariant operators from which it is possible to isolate the 
particular topological invariant used to fix the topological constraints. If the topological 
invariant is the Gauss linking number and the two topologically linked rings are represented 
as continuous curves embedded in the space and parametrized by their arc-lengths, this task 



has been achieved in 
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4l| . In the present case, the two rings are constructed out of a set 
of 2s open subtrajectories parametrized by the z coordinate and not by the arc-length. This 
parametrization is very peculiar because it identifies the parameter specifying the positions 



of the monomers with one coordinate of the space in which the monomers are 
For all the above reasons, the passage to topological field theories explained in 40 



uctuating. 



4l| cannot 



be straightforwardly applied to the present situation and has required a separate derivation. 
As a result of this derivation, we have been able to show that the path integral expressing 
the probability function of a 2s— plat formed by two polymer rings entangled together is 
equivalent to the correlation function of a gas of 2si particles of type 1 and 2s2 particles 
of type 2, where Si + S2 = s. The interactions between these particles are mediated by the 
vector fields of an abelian BF model. This is a topological gauge field theory that has been 
discussed in 36|,|42|. The particles are also subjected to short-range interactions whose ori- 



gin is the following. The 2s directed paths compo sing the two-polymer link are treated here 



like paths of directed polymers in random media 
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44| . which are subjected to quenched 



random potentials. After integrating over the random noise according to the prescriptions 



of Ref. J43|, in the polymer action of the 2s directed polymers appear potentials describing 
short-range forces acting on the monomers. 

The final passage to field theory is performed using the analog in statistical mechanics 
of the second quantization process. To this purpose, we generalized the method used by 
de Gennes and coworkers J45| to achieve the field theory formulation in the case of a single 
polymer chain subjected to short-range interactions to the case of a set of 2s different 
polymers. In the "second quantized" version of the statistical mechanics of the 2s— plat, the 
scalar fields create and destroy monomers in different positions of the space. The square 
module of such fields may be related to the monomer density at a certain point. The BF 
fields take into account the interactions necessary to keep the system in its initial topological 
configuration. The abelian BF model has been quantized in the Coulomb gauge, because in 
this gauge the analogy with anyon field theories becomes particularly explicit. The obtained 



field theory is a multicomponent model of anyons such those described for instance in j46|. 
This kind of theories exhibits the phenomenon of superconductivity. The only difference 
in our case is that the scalar fields containing the creation and annihilation operators for 
particles of type 1 and 2 are organized in replica multiplets, where at the end the limit of 
zero replicas should be taken. 

This paper is organized as follows. First of all, to map the partition function of two 
linked polymer rings into that of anyons, it is necessary to split their trajectories into a 
set of 2s subtrajectories, which in the anyon model describe the evolution in time of the 
quasiparticles. The splitting procedure and the definition of a time variable that is able 
to parametrize the 2s subtrajectories is carefully described in Section [TTl A proof that it is 
possible to isolate from the amplitudes of the BF model the Gauss linking number also after 
splitting the trajectories and changing their parametrization, is presented in Section [XTTl The 
fact that after quantizing the BF model in the Coulomb gauge it is still possible to recover 
the Gauss linking number from the amplitudes of the holonomies is shown in the particular 
case of a 4— plat in Appendix B. In Section [IV] the partition function of two linked polymers 
subjected to topological constraints imposed with the help of the Gauss linking invariant is 
transformed into a theory of 2s directed polymers interacting with the magnetic-like fields of 



the BF model. Contrarily to Ref. |28|, we treat the 2s subtrajectories as trajectories of real 
directed polymers. This requires the introduction of random potentials which complicates 
somewhat the passage to the anyon field theory performed in Section |Vl In the anyon 
formulation, the densities of monomers associated to the two original polymer rings can 
be regarded as the densities of anyon quasiparticles of type 1 and 2 interacting together. 
Thanks to a Bogomol'nyi transformation, the interactions may be split into a self-dual part 
and a part containing only short-range interactions. Remarkably, the latter interactions 
persists even if the short-range interactions coming from the random media are switched off. 
This is an effect of the presence of the topological constraint. In Section |VT] it is reviewed 



for completeness the case of a 4— plat studied in [28|. The static configurations of the anyon 
densities that minimize the Hamiltonian are computed. It is shown that the anyon model 
admits static self-dual points. The nature of the density configurations corresponding to 
these self-dual points is analyzed in Section IVIII We prove that the solutions of the classical 
equations of motion that minimize the static Hamiltonian are self-dual configurations, whose 
exact form can be obtained after solving a sinh-Gordon equation. Finally, our conclusions 
are drawn in Section IVIIII 
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FIG. 1. Representation of a treefoil 

knot in terms as a two-dimensional FIG. 2. The figure shows one of the 

diagram. crossings which are present in the di- 

agram of the treefoil knot of Fig. [1] 

II. POLYMERS AS 2s-PLATS 

Let's consider two closed loops Fi and r2 of lengths Li and L^ respectively in a three 
dimensional space with coordinates (r,^). The vectors r = [x^y) span the two dimensional 
space M^. z will play later on the role of time. The two loops will be labeled by using a 
indices the first letters of the latin alphabet: a,b, . . . = 1,2. We will assume that Fi and 
F2 form a 2s— plat. For convenience, we briefly review what is a 2s— plat. First of all, we 
recall that a single closed trajectory is from the mathematical point of view a knot, while 
a system of knots linked together forms a link. Knots and links may be represented after a 
projection onto a plane by diagrams like those of Fig. [1] and Ej in which the original three- 
dimensional structure is simulated by a system of crossings, see Fig. |2l Each crossing is 
composed by three arcs, one overpass and two underpasses. One may also realize that the 
treefoil diagram in Fig. [1] is characterized by two minima and two maxima. Two dimensional 
diagrams of this kind, deformed in such a way that the number 2s of minima and maxima 
is the smallest possible and the maxima and minima are aligned at the same heights ZMax 
and ZMin respectively, are called in knot theory 2s— plats u. The height of a 2s— plat is 
measured here with respect to the z axis. 2s— plats are used to classify knots and links by 
^ Actually, to be rigorous one should still require that neither maxima nor minima occur at the crossing 
points. 







FIG. 3. A link formed by two polymers Pi and P2 



dividing them into classes characterized by the same value of s. The concept of 2s— plats 
arises naturally in biochemistry, see e. g. 11|. In the present case, with some abuse of 
language, we will call 2s— plats also the two dimensional diagrams in which maxima and 
minima are not aligned, like for instance in Fig. |3l Let us denote with the symbols Tq/^, 
/q = 0, . . . , 2sa — 1, the heights of the maxima and minima of each trajectory Fq, for a = 1, 2, 
Arbitrarily, we choose Ti^q and T2,o to be the heights of the lowest minima on the trajectories 
Fi and F2 respectively. Starting from Tafl-, we select the orientation of Fq in such a way that, 
proceeding along the trajectory according to that orientation, we will encounter in the order 
the points Ta^i, Ta,2i ■ ■ ■ 1 'Ta,2sa- Clearly, Tq^i is a point of maximum, Ta^2 one of minimum and 
so on. Moreover, we should put for consistency: 



Ta,2sa = '^afl 



(1) 



The introduction of this double notation for the same height r^.o will be useful in the following 
in order to write formulas in a more compact form. In the following the 2s— plats Fi and F2 
will be decomposed into a set of directed trajectories Taj^-, a = 1,2 and /„ = 0, . . . , 2sa — 1, 
whose ends are made to coincide in such a way that they form the topological configuration 
of two linked rings. An example when s = 3 is presented in Fig. HI In the general case the 
set of points belonging to F^,,/^ can be described by the formula: 



'- a,Ia ~ \ ^a,Ia\^a,Ia 



a =1,2; Ia = l,. 


. . , 2s, 


Ta,Ia-l < ZaJa < TaJa 


la odd 


TaJa < ^a,Ia < '^a,Ia-l 


/„ even 



(2) 
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FIG. 4. Sectioning procedure for a 2s-plat Fa with s = 3. 



where the additional conditions: 



raA(V/J = ra,/,+l(V/J Ia=l,...,2Sa-l 



(3) 
(4) 



which connect together the subtrajectories Taj^ so that the loop F^ is reconstructed, are 
understood. In Eq. ([2]) raj^{zaj^) represents the projection of the trajectory Taj^ onto the 
plane x,y transverse to the longitudinal direction Zaj^. Let us note that we are using the 
same indexes la to label the trajectories F^/^ and the points Taj^. However, in the first case 
Jq = 1, . . . , 2sa, while in the second case we have chosen /„ = 0, . . . , 2sa — 1- In the case of 
the variables Zajjs, the range of /„ is the same as that of the F^^/^'s. 

We notice that the variables Zajjs are always growing and do not take automatically 
into account the fact that the whole chain is continuous and has a given orientation. Better 
variables, both with respect to the continuity and orientation, are the following ta,ijs: 



t 



a, I a 



taja = ^a,ia when la is odd 

-{ZaJ, - Taj J + TaJ,-l whcU /« is CVCU 



(5) 
(6) 



Assuming for instance that la is odd, for two consecutive trajectories Taj^ and Taj^^-^^, we 
have that: 



Ta,Ia-l < ia,Ia+ < V-fa 



(7) 
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and 



ra,Ia > ta,Ia+l > Ta,Ia + l (8) 

According to the above conventions, trajectories labeled by odd /^'s are oriented from a 
point of minimum to a point of maximum, while trajectories with even values of /„ go from 
a point of maximum to a point of minimum. In the new coordinate taj^, the trajectory F^ /^ 
becomes parametrized as follows: 



^a,Ia 



T^a,Iaita,Ia) 



a =1,2; /a = l,...,2sa 

Ta,Ia-l < ta,Ia < 'TaJa h odd 



(9) 



where the boundary conditions ([2]) and ^ are understood. 

III. THE GAUSS LINKING NUMBER AND THE ABELIAN BE EIELD 

THEORY 



To express the topological properties of the system of two linked loops Fi and F2, we use 
as a topological invariant the Gaussian linking number: 



x(ri,F2 



Ati^^"" 



dxi{ai) (p dx2{(J2 



Xi{oi) -X2{a2)y 



3^2(0-2 



(10) 



where the x^(cra)'s, a = 1, 2, are closed curves representing the loops Fi and F2 in the three 
dimensional space. The variables cJi and (72 used to parametrize Fi and F2 represent the 
respective arc-lengths of the two loops. They are defined in such a way that < cTq < L^. 
In the following, the trajectories of the two loops will be topologically constrained by the 
condition: 

m = x(ri,F2) (11) 

m being a given integer. The above constraints is imposed by inserting the Dirac delta 
function 5{m = x(Fi,F2)) in the partition function of the 2s— plat, where the statistical 
sum over all conformations of Fi and F2 is performed. Of course, the analytical treatment 
of such a delta function in a path integral is difficult. Some simplification is obtained by 
passing to the Fourier representation: 

dX 



5(m-x(Fi,F2))= f 



,-iA(m-x(ri,r2)) 



^2n 



(12) 
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However, even in the Fourier representation, the difficulty of having to deal with the Gauss 
linking number in the exponent appearing in the right hand side of Eq. (fT2|) remains. For- 
mally, this topological invariant introduces a term that resembles the potential of a two-body 
interaction. However, this potential is both nonlocal and nonpolynomial. It is for that reason 
that the treatment of the Gauss linking number in any microscopical model of topologically 
entangled polymers is usually very complicated. The best strategy do deal with it so far 
consists in rewriting the delta function 6{m — x(ri,r2)) as a correlation function of the 



holonomies of a local field theory, namely the so-called abelian BF-model 

/•+00 
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47|: 



(13) 



where 






(14) 



In the above equation we have put x = (x, t) to be dummy integration variables in the three 
dimensional space. Moreover, Sbf[B,C] denotes the action of the abelian BF-model: 



Sbf[B, C] 



An 



d'xB^,{x)d^Cp{x)e^"'f 



(15) 



e^'^^, fijVjP = 1,2,3, being the completely antisymmetric e— tensor density defined by the 
condition e^^'^ = 1. k is the coupling constant of the BF-model. Finally, the constants Ci 
and C2 are given by: 



ci = A 



C2 



K 

8^ 



(16) 



While there is some freedom in choosing Ci and C2, one unavoidable requirement in order 
that Eq. flT3l) will be satisfied is that one of these parameters should be linearly dependent 
on K. In this way, it is easy to check that k, may be completely eliminated from Eq. (fl^ by 
performing a rescaling of one of the two fields S^ and C^. This is an expected result, because 
K, did not appear in the left hand side of Eq. ( IT3l) . so that it cannot be a new parameter of 
the theory. By introducing the currents: 



J!:ix) 



dS:>^{aa)d^^\x - XaicTa)) a = 1,2 



(17) 



r„ 



Zbf{X) niay be rewritten in the more compact way: 
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In all the above discussion, the two trajectories Fi and r2 have been parametrized with 
the help of the arc-lengths cJi and a2. However, in the case of a 2s— plat, the loops F^ 
are realized as a set of open paths F^/^ connected together by the conditions ([3H1]). The 
subtrajectories F^j^'s are directed paths VajS'^a,ia) = (^a/a('^«.'fa)' -^a/aX^^.^a)) parametrized 
by the new variables taj^-, which are connected to the third spatial coordinates Zaj^ = x^ 
by the relations ([5]) and ([6]). Due to this difference of parametrization, the above method to 
express the Gauss linking number based on the BF- model, in particular Eq. (fT3l) . should be 
changed appropriately. Our starting point is the new partition function: 

where Sbf[B, C] coincides with the action (IT5l) but with the fields B, C renominated B, C 
and 

J^M = Cai2 r'^ </.(^a./J^^'H^ - ^a,/.(ta,7j) (20) 

7^ = 1 JT-a,Ia-l 

Let us note that the transformation Za^i^ — )• ta,/^ provided by Eqs. ^ and (j2]) leaves the 
BF action and the source terms unaffected, so that it does not change the form of the 
path integral ZbfW- Fot this reason, starting from Eq. (120|) . the directed paths F^/^ are 
parametrized with the variables taj^ instead of Zaj^ ■ 

To show that the partition function Zbf{^) of Eq- ( fT9|) coincides with the partition 
function Zbf{^) of Eq. ( ITSl) . we introduce new variables (Ja,ia ^^ follows. Let <^a,iai ^a = 
0, . . . , 2sa — 1 be the arc-length on F^ of the point of maximum or minimum located at the 
height Taj^- The Cajjs are defined in such a way that they span the intervals: 

<;a,Ia < (^a,Ia < 'ia,Ia + l (21) 

As a brief digression, even if it is not necessary for the present discussion, let us define the 
arc- length (j'^j^ of each trajectory Taj^- It is easy to show that cr'^^j^ is given by: 

^1,7. = (^a,ia - "iaja (22) 

As a matter of fact, a'^^ j^ ranges in the interval [0, <ia,ia+i ~ 'iaja] ^^^1 ?a,/a+i ~ ^a,/a is the total 
length of Fa /^ for a = 1,2 and la =, ■ ■ ■ , 2sa — 1. 

On each subtrajectory F^ /^ defined by Eq. (^, we can separately pass from the parameters 
ta,ia to the arc-length of Taj^ by a transformation of the kind taj^ = ta,iS(Ta,ia)- Putting 

KjA'^aM) = ^ajAta,lMa,lJ) (23) 
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we may rewrite the currents J^ (x) as follows: 

J!ii^) = Y.Ca r'^ dx^^^j^{aa,lJS^'\x - Xa,lMa,lJ) (24) 

1^ = 1 -^iaja-l 

It is easy to realize that the sum over all values of /^ in Eq. (1M|) is equivalent to a contour 
integration over the whole loop r^. As a consequence: 

J^{X) =cA dx^MS^'Hx - iaicTa)) ^ J„^(x) (25) 

JTa 

i. e. the currents J^ix) coincide with the currents J^{x) defined in Eq. (fT7|) . Therefore, 
in the partition function Zbf{^) of Eq. (IT9|) the currents Jl^{x)'s may be replaced by the 

J,^(x)'s: 

Zbf{X) = /pS^PC^e-^^^^I^'^l 

Due to the fact that the BF fields S^ and C^ are just dummy field configurations over which 
a path integration is performed, it is possible to rename them B^ and C^ respectively. In 
conclusion, starting from the partition function Zbf{^) of Eq. ( IT9l) . the partition function 
Zbf{^) appearing in Eq. ( IT3|) has been recovered. In other words, it has been shown that: 

Zbf{\) = Zbf{\) (27) 

and thus in the identity (fT3l) we can replace Zbf{,^) by Zbf{^)'- 

/ + 00 
dXe-'^'^ZBFiX) (28) 

-oo 

This is the desired final result. Thanks to Eq. (l28l) . it will be possible to transform the path 
integral over all conformations of the 2s— plat, which is complicated by the cumbersome 
presence of the dirac delta function containing the Gauss linking number, into a path integral 
over the trajectories of a system of 2s particles interacting with magnetic fields. 

In order to establish the analogy between polymers and anyons, which will the subject 
of the next Section, it will be convenient to quantize the BF model in the Coulomb gaugcl: 

diB' = diC' = (29) 



^ A similar approach like that proposed here can be found in [48^. In Ref. [48l| the plats are however static, 
they do not fluctuate, and the light-cone gauge has been used. 



14 
After the gauge choice (129|) . the action of the BF model ( !T5|) becomes: 

SbfMB, C] = ^ fd^x[Bse''d,C, + Cse*^^,^,] (30) 

with e^^ = e*-'^ being the two-dimensional completely antisymmetric tensor. The gauge fixing 
term vanishes in the pure Coulomb gauge where the conditions ( l29l) are strictly satisfied. 
Also the Faddeev- Popov term, which in principle should be present in Eq. (l30l) . may be 
neglected because the ghosts decouple from all other fields. Moreover, the requirement of 
transversality of ( l29|) in the spatial directions implies that the spatial components Bi and 
Ci of the BF fields may be expressed in terms of two scalar fields b and c via the Hodge 
decomposition: 

Bi = SijCPh Ci = Eij&'c (31) 

After performing the above substitutions of fields in the BF action of Eq. ( !30|) . we obtain: 

Sbf,cg = £ y d^x[BsAc + C3A6] (32) 

Let's compute now the propagator of the BF fields: 

G^,(x, t; y, t') = (S^(x, t), C,{y, t')) (33) 

From Eq. ( l30l) it turns out that only the following components of the propagator are different 
from zero: 

G3.(x, t; y, f) = ^-^^^e,,d^^ log |x - yp (34) 

G'i3(x, t; y, t') = -G3i(x, t; y, t') (35) 

The path integration over the scalar fields b and c in the partition function ZbfW is gaussian 
and could be in principle performed. A natural question that arise at this point is the 
interpretation of the topological constraint f lTT]) in the Coulomb gauge? As a matter of fact, 
the BF propagator in the Coulomb gauge breaks explicitly the invariance of the BF model 
under general three dimensional transformation. It seems thus hard to recover the form (fTOj) 
of the Gauss linking number in this gauge. Of course, an equivalent constraint should be 
obtained in the Coulomb gauge due to gauge invariance. In Appendix B it will be shown 
by a direct calculation in the case of a 4— plat that this is actually true. The computation 
of the expression of the equivalent of the Gauss linking number in the Coulomb gauge for a 
general 2s-plat is however technically complicated and will not be performed here. 
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IV. THE PARTITION FUNCTION OF A PLAT 



In order to write the partition function of a 2s— plat, we follow the strategy explained in 



the previous Section of dividing each trajectory r^ into 2sa open paths F^ /^, /„ 



2Sa. 



The statistical sum Z{X) of the system, that is performed over all possible configurations 
^a,ia{'t) of the subtrajectories Taj^ using path integral methods, is defined as follows: 

2 2Sa 



Z{m) = n n /boundary ^r.,,. (t,,,Je-(^/-^+^-)5 (m - x(ri, r^)) 

a=l Ia=i conditions 



(36) 



In the above equation the boundary conditions on the trajectories enforce the constraints 
([3]) and (jl]). The free part of the action Sfree is given by: 



2 2Sa 



s 



free 



'^aJa 



a=l Ia = l "^■^".^a-1 



rfra,/,(ta,/J 



f^a,/. 



(37) 



The parameters Qaj^^ with that Qaj^ > 0, are proportional to the inverse of the Kuhn lengths 
of the trajectories r^/^. They are also related to the total lengths of the trajectories Fq/^ as 
it is discussed in Appendix A. Let us note that Sfree is a positive definite functional despite 
the presence of the factors (—1)^""^. This can be easily proved by performing inside Sfree 
the transformations of Eqs. (jS]) and ([6]) from the tajjs to the Zaj^ variables: 



Sfree = E E / ^^Zajadaja 



dYajSZajA, 



dZaJa 



(38) 



It is now evident that Sfree is either positive or, if the ra,/^(-2a,i'a)'s are constants, equal to 
zero. 

Since we wish to stress the analogy with directed paths moving in a random media, 
we have also to introduce a contribution with short-range interactions coming from the 



integration over the random noises 



43J. This is the origin of the contribution Sev to the 



total action in Eq. (15^ . Sev is of the form: 

Sev 



= E E / ^^1.^ / rft2,j(-l)^+^-V(ri,,(ti./)) - r,,jit,,j))5{h,j - h,j) 

7=1 J=l "'■'"I,/-! •'f"2,J-l 

+^EE E r'" ^^'^'^^ r'" ^^'^.^.(-i)'"'''""'^(ra,/.(wj)-w.(iMj) 



a=l Ia = l Ja=l ""^aJa-l 



X d{ta,I, -ta,Ja] 



(39) 
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In the right hand side of the above equation, the first part describes the interactions between 
the monomers belonging to different loops, while the second part takes into account the 
interactions between the monomers of the same loop. If the random noises are gaussianly 
distributed, the two-body potential V^(r) is of the form: 

V^(r) ~ Vo5{r) (40) 

with Vq being a positive constant. The factors (— l)^"*""'"^ and (^—ly^+Ja-^ appearing in 
Eq. (139|) are necessary to make the interactions repulsive. This can be easily proved by 
passing to the variables Zaj^ using the transformations of Eqs. (|5]) and ([6]). 

As explained in the previous Section, the delta function 6{m — x(ri, ^2)) appearing in the 
partition function Z[m) of Eq. fl5^ may be simplified by introducing the BF- fields B^, C^ 
with action 

Sbf = ^ fd^xe^'^'B^d^Cp (41) 

After performing the Fourier transform of the delta function according to Eq. (IT2|) and 
exploiting Eq. (l28l) . the partition function Z{m) becomes: 

/+00 
d\e-'^'^Z{\) (42) 

-00 

where 

Z(A) = |D5,2)ae---[-'-] n ft /.o.naa. ^r,,.(4,Je- (43) 

1—1 /a=l conditions 

The polymer action S can be split as follows: 

S = Sfree + SeV + Stop (44) 

The expressions of Sfree and Sev are provided in Eqs. ( 137|) and ( l39l) respectively. Finally, 
using Eq. (fT9|) . it is possible to realize that the topological contribution Stop to the polymer 
action turns out to be: 

Stop = 1 f d'x [J^ix)B^{x) + J^ix)C^{x)] (45) 

where the currents Ja{x), a = 1,2, are given in Eq. (|20|) . 
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V. AN ANYON FIELD THEORY FORMULATION OF POLYMERIC 2s-PLATS 

The starting point in this Section is the polymer statistical sum Z{X) of Eq. (I43l) . This is 
formally equivalent to the partition function of a multicomponent system of anyon particles. 
To write this partition function in terms of fields, we have to perform an integration over 
all polymer trajectories Taj^itaj^), a = 1,2 and la = 1, . . . ,2sa- The passage to the field 
theoretical formulation is not just a formal step, it allows to describe the short range and 
topological interactions by means a local and polynomial action. Before the introduction of 
fields, these interactions are both nonlocal and nonpolynomial. The standard procedure in 
polymer physics to pass from polymer trajectories to monomer densities and thus to a field 
theory consists in introducing auxiliary fields. In the case of the topological interactions, we 
have already seen that the auxiliary fields are the BF fields B^{x) and C^{x). The short 
range interactions in Eq. (139|) require instead several scalar fields in order to be simplified. 
The minimal number of these fields is 2sa + 2sb + 2. A couple of fields 01 (x, t) and 02 (x, t) is 
needed for the interaction between monomers belonging to different loops. The interactions 
between monomers belonging to the same loop will be taken into account by the fields 
V9i,/^(x,t)'s and (/?2,Jt,(x, t)'s with la = l,...,2sa, Jb = l,...,2sb. 



The passages that lead to the final field theory are well known in the polymer litera- 



ture 



le p 
4q|, 



41 



49 



50|. After an integration over the auxiliary scalar fields 0i(x, t), 02(x, t) 
and (^i^/^(x, t), (/?2,Jt,(x, t), /„ = l,...,2sa, Jb = l,---,2sfo, whose details are explained in 
Appendix C, the expression of the polymer partition function Z{X) of Eq. ( l43l) becomes: 



Z{X) = hin fvB^VC^ II /"p^i,,I)^t,/^M(ri,/(n,/-i),ri,7_i)V^}}(ri,,(ri,7),ri,,) 

"2^0 -^ \_I=l -^ 

2S2 „ 

J 1 J 



where the ^^/^ and ^ajai (^ = 1? 2, 7^ = 1, . . . , 2sa, are complex rephca fields: 



^a,i. = «,., . . . , C/j K,i. = «/a' • • • > c:/j (47) 
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The action A in Eq. (H6!) contains the free part and the topological interactions: 



"^ij 



1=1 "'''"1,1-1 



/-I 



(fx^lj{x,t) 



d 



dt 4(71,7 



\i-i 



V - 2A(-l)'-^B(x, t) + tX{-iy-'Bsix, t) 



^M(x,t) 



2S2 



J=l >^T2,J-l 



d 



IK 



dt 4(72, J V 27r' 
The action AeV) given by 



l)^-^C(x,t) + 



ZTT 



\J-1 



Cs{x,t) 



^2„/(x,t) (48) 



^ 



EV 



^ 2 

/,/' = ! 


/ dt 

'T\,I-\ 


2S2 

^ 2 

j,.j'=i 


/•T2..J 

/ (it 

-^ T2,J-1 



(it'(5(t - t') / (/2x(-l 



n/+/'-2 



^l,7(x,t) ^l,/'(x,t' 



cit'(5(t - f) / (i2x(-l 



\J+J'-2 



*2,j(x,t) 



^..7'fx,t'^ 



2si 2s 



+ ^oEE i' ^^ r' ^^' /rf'x(-i)^+^-2|^i,,(x,t)n^2,j(x,t')i''^(i-o (49) 

7=1 J=l "'''"1,-f-l >'T2,J-l "^ 

is the analog of the action 5*^^ written in the language of second quantized fields and 
describes the short-range interactions. Looking at Eqs. ( H6|) - (H9|) . we see that we have suc- 
ceeded in our task, i. e. the original polymer partition function ( H3l) has been transformed 
in an anyon field theory. The action A is formally equivalent to the action of a multicom- 
ponent system of anyons subjected to the Coulomb interactions described by Aev- Similar 
systems have been discussed in connection with the fractional quantum Hall effect and high 
Tc superconductivity 46|]. The only differences are the boundaries of the integrations over 



the time, which in the present case depend on the heights of the points of maxima and 
minima of the two trajectories Fi, r2 and the fact that here the quasiparticles are bosons of 
spin ni or n2 considered in the limit Ua — ;• 0, a = 1,2. 

VI. SELF-DUALITY OF THE TWO-POLYMER PROBLEM 



In this section we restrict ourselves for simplicity to 4— plats. The partition function of 
a 4— plat formed by two linked polymers is obtained by putting si = S2 = 1 in the general 
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partition function of a 2s— plat given in Eq. (H6l) . Accordingly, the action A of Eq. (|38 
becomes in this particular case: 



n,i 
A= I dt I rf2x<j ^*^ 



9 



at 4(?i, 



-D2(-A,B) + aS3 



^1,1 



+ ^* 



<9 



dt 4(71, 



-D2(A,B)-iA53 



^1, 



T2,l 



+ / dt I rf^x^;,i' 

't2,0 



a 



dt 4^2,1 



D^ 



-^■°)-^-= 



^2,1 



+ ^ 



2,2 



« l.D^f^.c)-|c3 



^ 



dt 4c/2,2 V27r 
In writing the above equation, we have used the notations: 



2,2 



(50) 



D(±A,B) = V±zAB 



b(±—,c 

V 27r' 



K 



V±z— C 

27r 



(51) 



The short-range interaction term Aev oi Eq. ( 149|) simplifies in the case of a 4-plat as follows: 

2 /■min[Ti,i,T2,i] /• 

^W= 5Z ^0 / rft /d2x|^i,,(x,t)n^2,j(x,t)|2 

72„hf, ^„ +M2|,f, ^„ +M2 



Vn 



-E? 



2 



dt /d^x|^i,,(x,t)H^i,,,(x,t)|^ 



+ E Y r''rfi/^'x|^2,j(x,t)n^2,j'(x,t)p 



(52) 



Finally, the BF contribution iSbf defined in Eq. (HT]) remains unchanged. 

The next goal is to find the classical field configurations which minimize the the energy 
J-" of the two-polymer system. A sketchy derivation of these configurations can be found in 
Ref. |28j]. In the following, we will provide the details that were missing in 28|. The energy 
J-" is given by: 

jr = iS^F + A + Aev (53) 

where the expressions of Sbf, -^ and Aev are defined in Eqs. (141 p . (I50p and (152]) respec- 
tively. To simplify the task of its minimization, the short-range interactions will be neglected 
putting Vq = in Eq. (I52p . so that Aev = 0. This approximation is valid for instance for 
polymer solutions which are at the theta point. To proceed, we notice that the third compo- 
nents i?3 and Cs of the BF fields play the role of pure Lagrange multipliers. Thus, they can 
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be integrated out from the partition function (H6l) giving as a result the following constraints: 



B = 2(|^2,i|' - |^2,2r)^(r2,i - mt - r2,o) 



C 



AttX 



K 



(|M/ 



i,ir-|^i,2r)^(ri,i-t)^(t-ri,o) 



(54) 
(55) 



where B and C are the magnetic fields associated to the vector potentials Bi and Cj respec- 
tively: 

B = diB2 - 8281 = e'W.Bj (56) 

C = diC2 - d^Ci = e'^diCj (57) 

In Eqs. ([511) and ([55]), 9{t) denotes the Heaviside function 6{t) = if t < and 9{t) = 1 
if t > 0. We will look here only for static field configurations, i.e. those which satisfy the 
relations: 

|<1 = |€?.=0 (58) 

for all values of a = 1, 2 Jq = 1, 2 and 0"^ = 1, . . . , ria, where the n^s denote the numbers of 
replicas. To avoid problems with the presence of the Heaviside functions in the expression 
of the magnetic fields, we will assume that 



Tl,0 = ^2,0 = To 



and 



Tl,l = T2A = Ti 



(59) 



In this way, the parameter t, whose range is changing depending on which subtrajectory is 
parametrized, is always be defined in the interval [tqjTi] as a real time. At this point, the 
static energy J-'st may be written as follows: 



J^. 



St 



[n - To) 



+ [ri - To) 






1 



4<7i,i 
1 



D(-A,B)^i,i 



+ 



1 



4^1,2 



D(A,B)^i,2 



4^; 



2,1 



D(--,C)vl/2,i 



+ 



1 



4^2,2 



D(^,C)v|/2,2 



(60) 



The vector potentials B and C in the above equations are determined by the relations 
f l5l] - [571) . The analogy with the anyon problem suggests the application of the Bogomol'nyi 



identities 



51] . For a single theory of complex scalar fields i/j* , ip minimally coupled to an 



abelian gauge field a, these identities look as follows: 

|D(7, a)^|2 = |D±(7, a)^|2 ^ 7b|^|2 ± e'^dijk 
where D(7, a) = V — i^Si is the covariant derivative and 



(61) 



D±(7,a) = Di(7,a)±«D2(7,a) 



(62) 
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Here .0,(7, a), i = 1,2, denotes the components of D(7, a), while 



b = dia2 - d2ai 



(63) 



is the magnetic field. Finally 



jk = — [^*/5fc(7, a)?/^ - ^Dk{-f, a.)^*] 



(64) 



is the current related to the abelian gauge group of symmetry. Let us notice that the term 
in Eq. (16T]) containing j^ is a total derivative, so that it can be omitted in our case, in which 
the space has no boundaries. 

Coming back to the problem of minimizing the static free energy of Eq. f l60ll . we can now 
apply the Bogomol'nyi identities for all 2ni + 2n2 replica fields. Actually, Eq. ( 16T1) defines 
two different identities, depending on the choice of sign. This fact may be used to simplify 
the calculations. In particular, we will choose the + sign when the scalar fields are coupled 
to B and the — sign when the scalar fields are coupled to C. As a result we obtain: 



J^st — J^i + J^2 



(65) 



with 



^1 



+ 



[Tl - To 



D 



) d^< 



^9 



1,1 



4^2,1 



K 



-^,C)^,, 



D+(-A,B)^i 
1 



2 1 
+ 



Ag 



1,2 



D+(A,B)^i,2 



+ 



4c/: 



2,2 



D_(— ,C)^2,2 



(66) 



^2 = in - To) / d^ 



A 



4^ 



-B 



1,1 



^1, 



4^ 



-B 



1,2 



^1, 



K 



i'^92,1 



-C 



^2,1 



+ 



K 



8^92;. 



-c 



^ 



2,2 



(67) 



If we substitute in Eq. (167|) the expressions of the magnetic fields B and C in terms of the 
replica fields given by Eqs. (1M|) - (155]) . it is easy to realize that J-2 becomes of the form: 



^2 = -(n -To) / rf^x 



1 
91,1 



ll^2,l| ~ l^2,2| ) |^l,l| ll^2,l| — I '±'2,2 1 ) |^1,2| 



1 
9l,2 



f |^l,l|' - 1^1,2!') |^2,l|' + ( l^l.ll' - 1*1,2 

52,1 V / g2,2 



(68) 
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It turns out from the above equation that the presence of the topological constraints induces 
changes in the energy of two linked polymers which consists in the appearance of short-range 
interactions with coupling constants proportional to 

^ (n-ro) (69) 



9a,Ia 



These interactions clearly interfere with the short-range interactions given in Eq. ( 152|) . which 
have potentials of the same structure, characterized by fourth-order powers of the fields, but 
have different coupling constants. In particular, in Eq. (l52l) the coupling constant Vq is 
always positive, while the coupling constants in Eq. (169|) can be either positive or negative. 
This shows that the topological constraints have nontrivial effects on the short-term inter- 
actions acting on the monomers. These effects have been already observed in experiments. 



see for example Ref. 



14| . Analytically, the influence of the topological constraints has been 



52 



54| . Thanks to the analogy be- 



quantitatively described using various approximations 
tween anyons and 2s— plats established here, we have been able to derive Eq. (168|) . which 
represents a direct confirmation at a nonperturbative level of the appearance of interactions 
associated to topological constraints. 

Let us now go back to the expression of the static energy J^^t of Eq. (165|) . Looking at 
the form of its components Ti and J^2 of Eqs. (166|) and ( 168|) . it is possible to conclude that 
J^st is formally equivalent to the Hamiltonian of a set of complex scalar fields coupled to 
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the BF fields B^ and C^. This kind of theory is known to have self-dual solutions 
In general, the search of self-dual solutions is not a simple task, because of the non-linear 
character of the classical equations of motion. Up to now, this problem has been solved in 
general only using numerical methods. Despite these difficulties, however, it is still possible 
to investigate the self-dual point analytically by restricting ourselves to the region of the 
space of physical parameters in which the attractive and repulsive forces appearing in ^-"2 
counterbalance themselves. In the present context, the self-duality is achieved when the 
following conditions are satisfied: 

9i,i = 9i,2 = 92,1 = 92,2 = 9 (70) 

If the equalities in Eq. (170|) are valid, in fact, the potentials in the right hand side of Eq. (168!) 
vanish identically. As a consequence, the energy (|65|) becomes self-dual, i.e. it can be written 
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as a sum of self-dual contribution: 

in - To 



^,. 



4<7 
4<7 



d'x 



d^x 



D 
D 



-A,B)^. 



+ 



K 



ZTV 



D+(A,B)^i,2 



D_(— ,C)^2,2 



(71) 



In anyon field theories the self-duality condition (170]) has a very physical meaning, see 
for example J5l|, but its interpretation in the case of the 2s— plat is much more difficult. 
Certainly the self-duality condition fITOl) is related both to the length and rigidity of the 
polymer trajectories. Indeed, the parameters gaj^ can be identified with the inverse of the 
Kuhn lengths of the subtrajectories Taj^ and thus determine their rigidity. Moreover, in 
Appendix A it is shown how the lengths of the subtrajectories depend on the Qaja^i ^^^ 
Eq. flAlSp . Therefore, it is clear that the relations (1701) are also imposing conditions on the 
lengths of the trajectories Fi i,ri 2,r2,i and r2,2, which must have in the average the same 
lengths in order to attain the self-dual point. 

In the next Section we will derive some explicit self-dual configurations which minimize 
the free energy J^st of Eq- (jTTl) . 



VII. SELF-DUAL SOLUTIONS OF THE TWO-POLYMER PROBLEM 



The task of this Section is to find classical solutions of the equations of motion which 
minimize the energy J^st of Eq. ( 1711) . The classical equations of motion read as follows: 

D+(-A,BX;i = (72) 

D+(A, 3)^^^,2 = (73) 

D-(-|^,c)v^2l = (74) 

D_(-|-,C)V^21 = (75) 

cri,cr2 being replica indexes. To Eqs. (172]475|) one should add the constraints (15^ and (155|) : 



e*^a,5, = 2 1^2,1 



e'^d^C, 



47rA 



|^2,2|'j^(Tl-t)0(t-To) 

|2 |,f, |2\ 



tv 



(|^i,ir-|^i,2r)^(ri-t)0(t-ro) 



(76) 

(77) 



To avoid analytical complications due to the presence of the Heaviside theta functions, we 
have assumed as in Eq. (159|) that tio = T20 = tq and ri^i = T21 = ti. Moreover, in the 
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following we will restrict ourselves to the replica symmetric solutions of Eqs. ( !72]475|) and 
(1761 1771) by putting: 



tjj^\ = ip^j for l<ai<ni I = 1,2 



^: 



2,J 



'4'2,J for 1 < (72 < 712 J = 1,2 



(78) 



In this way, the explicit form of the equations of motion ( 172H751) and of the constraints 
(17611771) looks as follows: 



[^1 - iXBi +i{d2- tXB2)] V'l,! = 
[di + iXBi + i{d2 + iXBi)] ^1.2 = 

V^2,i = 



^1 + — Ci -i[d2 + —C2 



^: 



2,2 







(79) 
(80) 
(81) 

(82) 



e^^a,5,- = 2ni(|^2,i|'- 1^2,2!') 
4?7,2vrA 



e'^diCj 



K 



,ir-i^i.2n 



At this point we pass to polar coordinates by performing the transformations: 



(83) 
(84) 



^^ 






(85) 



After the above change of variables in Eqs. (17911841) . we obtain by separating the real and 
imaginary parts: 



^iCJi^i - XBi + -a2logpi,i = 
-920;!, 1 + XB2 + -9ilogpi,i = 

(9iWi,2 + XBi + -92logpi,2 = 
-d2UJi,2 - XB2 + -di logpi_2 = 

K 1 

dlUJ2,l - 7— Ci - -92logp2,l = 

Ztt Z 

K 1 

d2UJ2,i - 7—^2 + -dx\o<gp2.\ = 

ZTT Z 

K 1 

dlUJ2,2 + Tj-C*! - 7:^2 logp2.2 = 

ZTT Z 

K 1 

d2UJ2,2 + TrC2 + -di logp2,2 = 

Ztt Z 



(86) 
(87) 
(88) 
(89) 
(90) 
(91) 
(92) 
(93) 
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e'WiBj = 2ni (p2,i - P2,2) (94) 

e'^d.C, = ^^^ (pi,i - P1.2) (95) 

K 

To solve equations ([HSHS3]) with respects to the unknowns Uaj^ and Pa,ia> ^^ proceed as 
follows. First of all, we isolate from Eq. (15^ and Eq. f lSS]) the same quantity XBi. By 
requiring that the expressions of A-Bi provided by Eqs. (1861) and (1881) are equal, we obtain 
the consistency condition: 

5iwi,i + 2^2 log pi,i = -(9ia;i_2 - -82 log pi,2 (96) 

A possible solution of Eq. ( p6|) is: 

wi,i = -t^i,2 and pi,i = ^- (97) 

Pl,2 

where Ai is a constant factor. As well, we could require that the two different expressions 
of the quantity XB2 obtained from Eqs. (187|) and ( 189|) are equal. However, in this way one 
obtains once again the condition ( l96l) . which can be solved by applying the ansatz (l97l) . In 
a similar way, it is possible to extract from equations ( l90] - [93l) the conditions: 

A2 
uj2,i = -UJ2,2 and p2,i = (98) 

P2,2 

with A2 being a constant. 

Thanks to (I97p and (198 p . the number of unknowns to be computed is reduced. For in- 
stance, if we know the expressions of cji^i,c<;2,i,pi,i and p2,i, the classical field configurations 
^1,2,^2,2, p{^,'^) and p2,2 can be derived using Eqs. (|97|) and (p8|) . As a consequence, the 
system of equations (l86l - [95|) reduces to: 

XBi = diui^i + -82 logpi,i (99) 

XB2 = 52a;i,i - -di logpi,i (100) 

K 1 

--Ci = diUJ2,i - -82 log P2,l (101) 

271 Z 

K 1 

--C2 = 52^2,1 + -di logp2,i (102) 

diB2 - d2B, = 2ni ( p2,i - — ) (103) 

H V P2,lJ 
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where we have used the fact that e^^diBj = diB2—d2Bi and e^WiCj = diC2—d2Ci. Eqs. fll04p 
contain only the unknowns cji,i,u;2,i,Pi,i and p2,i that have still to be determined. 

By subtracting term by term the two equations resulting from the derivation of Eqs. ( l99|l 
and f llOOp with respect to the variables x^ and x^ respectively, we obtain as an upshot the 
relation: 

A {diB2 - 92-Bi) = did2UJi^i - d2diUJi^i - -Alogpi,i (105) 

with A = 5^ + 5| being the two-dimensional Laplacian. 
Assuming that Wi i is a regular function satisfying the relation 

did2UJiA — d2diuji^i = (106) 

Eq. (llOSp becomes: 

X{d,B2-d2B,) = -iAlogpi,! (107) 

An analogous identity can be derived starting from Eqs. (llOip and (I102p : 

- idiC2 - d2Ci) = A log p2,i (108) 

The compatibility of (I107p and (llOSp with the constraints (I103p and (I104p respectively leads 
to the following conditions between pi i and p2,i: 

A log pi,i = AXm ( — - p2,i ) (109) 

A log p2,i = 4An2 ( pi,i - — ) (110) 

V Pi,i/ 

The fact that pi^i and p2,i appear in a symmetric way in Eqs. (llOOp and (IllOp . suggests the 

following ansatz: 

P2,l = -^ (111) 

Pl.l 

y4.3 being a constant. It is easy to check that with this ansatz Eqs. (I109p and (IllOp remain 

compatible provided: 

^ = -^ and ii = _!^ (112) 

^3 ni Ai rii 

We choose Ai to be the independent constant, while A2 and A^ are constrained by Eq. (I112p 

to be dependent on Ai. 

A2=(^yA^ A, = -^A, (113) 

\nij ni 
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We are now left only with the task of computing the explicit expression of pii. This may 
be obtained by solving the equation: 



[A, \ 

Pi,i 

VPi,i / 



A log pi,i = 4An2 -^ - pi,i (114) 

The other quantities p2,i, Pi, 2 and p2,2 can be derived using the relations fillip . ( !97|l and f l98|) 
respectively. Eq. (11141) may be cast in a more familiar form by putting: 77 = In ( -^^ j . After 
this substitution, Eq. (11141) becomes the Euclidean sinh-Gordon equation with respect to i]: 

A?7 = 8An2V Ai sinhr/ (115) 

Next, it is possible to determine the magnetic fields B and C from Eqs. (I103p and (I104p . 
In the Coulomb gauge, in fact, the two dimensional vector potentials B and C can be 
represented using two scalar fields b and c as follows (see also Eq. fl3T|)): 

B = (-^26, d^b) C = {-d2C, d^c) (116) 

Performing the above substitutions in Eqs. (11030 and (I104p . it turns out that b and c satisfy 

the relations: 

A6 = 2n2(pi,i - — ) (117) 

Pi,i 

4n27rA Ai 
Ac= (pi.i ) (118) 

1^ ' Pi,i 

The solution of Eqs. (I117p and (IllSp can be easily derived with the help of the method of 

the Green functions once the expression of pi^i is known. Finally, the phases a;i_i, a;i,2j ^2,1 

and a;2,2 are computed using Eqs. (p9|) -f lT02|) . In fact, remembering that we assumed that 

oji^i = —uji^2 and 002,1 = ^2,2 in (jHTj) and ( !98|) respectively, we have only to determine oJi^i 

and 0^2,1 . By deriving Eq. (p9l) with respect to z^ and Eq. (llOOp with respect to x^, we 

obtain: 

XdiBi = 9iWi,i + -8182 logpi,i 

Xd2B2 = dluji^x - -8281 logpi,i (119) 

On the other side, by adding term by term the above two equations and using the fact that 
in the Coulomb gauge the magnetic field B is completely transverse, it is possible to show 
that: 

Awi,i = (120) 
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Proceeding in a similar way with Eq. (llOip and f ll02p it is possible to derive also the relation 
satisfied by cj2,i: 

Awa,! = (121) 

VIII. CONCLUSIONS 

In this work a system of two polymers forming a nontrivial link has been considered. The 
topological properties of the link have been described by using the Gauss linking invariant. 
This is a weak topological invariant, but when applied to a 2s— plat configuration, which 
cannot be destroyed because the 2s points of maxima and minima are kept fixed, its capabil- 
ities to distinguish the changes of topology are greatly enhanced. The reason is the synergy 
between the constraint imposed by the Gauss linking number and those imposed by the fact 
that the polymer system cannot escape the set of conformations allowed in a 2s— plat. We 
have also seen in Appendix B what is the meaning of these constraints from the point of 
view of the 2s— plat. Basically, the sum of the winding numbers of all pairs of the subtra- 
jectories T^j^ are constrained to be equal to some integer multiple of 27r. Moreover, since 
the endpoints of the trajectories are fixed, also the winding number between two different 
trajectories is fixed up to multiples of 27r. Allowed are only the topology changes such that 
an amount of the winding angle of two subtrajectories is transferred in units of 2it to the 
winding angle of another couple of subtrajectories. This result paves the way to a treatment 
of polymer knots or links constructed from tangles. Polymers of this kind are relevant in 
biochemistry because nontrivial knot configurations appearing as a major pattern in DNA 
rings are mostly in the form of tangles |ll| . 

A crucial point of the connection shown in this paper between 2s— plats and anyons is the 
possibility of elimininating the cumbersome topological constraint (ITT]) from the partition 
function Z{X) of a 2s— plat by introducing BF fields. Indeed, the delta function fixing the 
constraint (ITT]) can be represented using the Fourier transform of the amplitude of gauge 
invariant and metric independent observables of an abelian BF model. This has been proved 
in Eq. ( 128|) . The proof of this relation is not trivial because the 2s subtrajectories in which 
the original 2s— plat has been split are open and parametrized by a special variable, the 
parameter t, which is proportional to one of the spatial components of the subtrajectories 
themselves. 
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Thanks to the identity (l28l) it has been possible to interpret the problem of the statis- 
tical mechanics of a 2s— plat as that of a two-component anyon gas with 2si particles of 
kind 1 and 2s2 particles of kind 2 interacting via short range potentials, see Eqs. (14214451) . 
The trajectories of the quasi-particles correspond in the polymer analog to the 2s directed 
trajectories Taj^-, a = 1,2, la = l,...,2sa, which are traversed by the fictitious currents 
(pOj) . The Gauss linking number can be interpreted as the circulation of the magnetic field 
generated by the current traversing the loop Fi with respect to the closed contour formed 
by the loop r2 
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55|. 



The system of quasiparticles with partition function Eq. (143|) has been further mapped 

into a two-component anyon field theory, whose final form is displayed in Eqs. (H6l) and (HHl - 

Q^. Similar field theories have been proposed in the past to explain the supeconductivity of 



high temperature superconductors without breaking the P and T invariance, see 



31|. The 



analogy between directed polymers and vortex lines has been studied in connections with 



high Tc superconductors in [56|. As in superconductors of type II, also in the present case 
attractive and repulsive forces appear, which vanish at some self-dual point of the theory. 
What is remarkable here, is that these interactions do not need the introduction of any 
potential and are purely related to the topological constraints imposed on the trajectories 
of the original polymer system. Indeed, they remain even if the short-range interactions are 
switched off as shown in Eq. (168|) . From the condition (170|) . which determines the existence 
of the self-dual point or not, it is possible to predict that 2s— plats consisting of homogeneous 
polymers should have a profoundly different behavior than their counterparts built out of 
block copolymers. As a matter of fact, we have seen that the physical characteristics of the 
2s directed polymers into which the 2s— plat has been split are described by the constants 
Qaja- ^^ particular, the rigidity of the trajectories may be specified by choosing the g'^j s 
appropriately. Clearly, from Eq. (170|) it turns out that the self-dual point is attained only 
if these constants are all equal, implying that either the polymer rings are homopolymers 
or their subtrajectories F^ /^ contain monomers of different types but, after averaging over 
the distance of many monomer sizes, they have identical physical properties. We have also 
derived the equations of motion that minimize the action the anyon field theory in the case 
of a 4— plat. These equation describe the self-dual point of the two-component anyon gas. 
After many simplifications, the relevant degrees of freedom can be derived by solving the 
sinh-Gordon equation ( I115p and the Laplace equations (I120H121]) . 
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From the polymer point of view, the physical meaning of the self-duality is unclear, 
because here only static conformations have been considered. In principle, these static 
solutions could become physically relevant in the case of a very long 4— plat in which the 
monomer concentration does not depend on the height z. What is however more important, 
is that Eq. flllSp . which determines the static density of monomers of type 1, is a sinh- 



Gordon equation identical to that obtained in [57| for the static vortices of a relativistic 
abelian Higgs model on a special type of Riemann surfaces. This analogy between field 
theories on Riemann surfaces and polymers, together with the connections between linked 
polymers and multicomponent anyon systems, that are related both to topological quantum 
computing and to high-Tc superconductors, should be further explored. It is true that for 
topological computing nonabelian anyon systems are necessary, while our discussion has been 
limited to the abelian case. However, this limitation is only apparent. In principle, instead 
of the Gauss linking invariant , we could have used much more refined topological invariants 



that would have led to nonabelian anyon field theories 
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. Up to now, however, nobody 



has succeeded to formulate completely a nonabelian theory of topological entanglement for 
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polymer systems based on such topological invariants, apart from a few exceptions 
Also the possibility of studying the statistical mechanics of knots constructed from tangles 
should be investigated, because up to now there is no analytical model which is able to 
describe the statistical properties of knots. 
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Appendix A: The length L of a directed polymer as a function of the height 

In this Appendix we consider the partition function 

Z = IVr{z)e-^ (Al) 

where S is the action of the free open polymer, whose trajectory F is parametrized by means 
of the height z defined in some interval [tq, n]: 



S = 9 



n 



dz 



To 



dr 



dz 



(A2) 



We want now to determine how the total length of the curve F depends on the constant 
parameter g. To understand what we mean by that, let us consider the standard case of 
an ideal chain whose trajectory is parametrized with the help of the arc-length a. We 
denote with a the average statistical length (Kuhn length) of the A^ segments composing 
the polymer. In the limit of large A^ and small a such that the product Na is constant, the 
total length L of the polymer satisfies the relation 



Na 



(A3) 



We wish to obtain a similar identity connecting L with A^ and g in the present situation, 
which is somewhat different. To this purpose, we first dicretize the interval of integration 
[to, ti] splitting it into A^ small segments of length: 

Tl -To 



Az 



N 



As a consequence, we may approximate the action as follows: 



N 



S ^ g^ 



■w=l 



Ar, 



Az 



Az 



where the symbol Ar^ means 



and 



(5r„ 



'■w+l 



(A4) 



(A5) 



(A6) 



r^ = r(ro + wAz) (A7) 

The discretized partition function becomes thus the partition function of a random chain 
composed by A^ segments: 

(A8) 



2m 



disc 



w=l 



^ lA |2 

r,„e "=i 



Az 
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Using simple trigonometric arguments it is easy to realize that the length of each segment 
is: 

AL = v/|Ar^|2 + A22 (A9) 

This is of course an average length, dictated by the fact that, from Eq. (1A8I) . the values of 
|Ar^„| should be gaussianly distributed around the point: 

Az 



|Ar„ 



9 



(AlO) 



In the limit Az — )■ 0, the distribution of length of Ar^„ becomes the Dirac (5-function: 



hm lj9-,\Ar^\yA. ^ ^ 1 1^ I _ Af 



Az^o 2y Az \ ' \l 9 

If A^ is large enough, we can therefore conclude that the total length of the chain F is 



(All) 



'Az 
L^ NAL = Nxl — + Az2 

9 



Since A^A^; = Ti — Tq, we get: 



L' = \n-To\' + 



N{n - To) 



(A12) 



(A13) 



In the limit A^ — )■ oo, while keeping the ratio — finite, Eq. ( 1A13I) becomes the desired relation 



between the length of T and g which replaces Eq. (\A3\i . 



Appendix B: The expression of the Gauss linking invariant in the Coulomb gauge. 



To fix the ideas, we will study here the particular case of a 4— plat. In the partition 
function (H3|) we isolate only the terms in which the BF fields appear, because the other 
contributions are not connected to topological constraints and thus are not relevant. As a 
consequence, we have just to compute the following partition function: 



Z, 



C'S,CG 



(A) 



VB^VC^.e"'''^^'^^ 



-Stot 



(Bl) 



where the BF action in the Coulomb gauge Sbf,cg has been already defined in Eq. (I3UJ) and 
Stop has been given in Eq. ( l45l) . In the case of a 4— plat. Stop becomes: 



s. 



top 



i\ / dt 

in r^'i , 



— ^5^(ri,i(t),t) j^5^(ri,2(t),t) 

dx2,iit) 



dt 



dt 
C,{v,^,{t),t)--^C,{v,,,{t),t) 



(B2) 
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where we recall that x^j{t) = {raj{t),t), a = 1,2, I = 1,2. Using the Chern-Simons 
propagator of Eqs. (1M|) - (135|) . it is easy to evaluate the path integral over the gauge fields in 
Eq. (IB II) . The result, after two simple Gaussian integrations, is: 



Zbf,cgW = ^^P ^ ^ X^ ( 
^ i,j=i 



'iy^'-%. 



dixUt) - xUt)) 



TO 






(B3) 



In the above equation we have put for simplicity: 



To = max[Ti,o,T2,o] 

Ti = min[ri,i,r2,i] (B4) 

For instance, if the polymer configurations are as in Fig. [5l we have that Tq = ri_o and 
Ti = ^2,1- Moreover, we remember that in our notation Tajit) = (x^ j-(t), x^ /(t)). Apparently, 




FIG. 5. Example of configuration of a 4— plat. 

the elements of the trajectories Fi and F2 which lie below Tq and above Ti do not take the 
part in the topological interactions. Thus is due to the presence of the Dirac (5-function 
6{t — t') in the components of the Chern-Simons propagator (1Ml)-f l35l) . However, we will 
see later that also the contributions of these missing parts are present in the expression of 
ZqsW- In order to proceed, we notice that the exponent of the right hand side of Eq. (IB3P 
consists in a sum of integrals over the time t of the kind: 

»T1 



-Di,/;2,j(ti) - Dij.2,j{to) = Si, 



TO 



d {x\,{t) - xlAt)) '-^ ^ (B5) 
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The above integrals can be computed exactly. It is in fact well known that the function 

Dij-2,j{'t) is the winding angle of the vector ri /(t) — r2^j(t) at time t: 

( x\j{t) -x\j{t)\ 
D,,,,At) = arctan [^^r^^^^^^^j W 

Thus, the quantity -Di,/;2,j(''"i) — Di,I;2,j{to) is a difference of winding angles which measures 
how many times the trajectory Tij turns around the trajectory r2,j in the slice of time 
T~o ^t < Ti. At this point, without any loss of generality, we suppose that the configurations 
of the curves Fi and r2 is such that the maxima and minima Taj are ordered as follows: 

T2,0 < 7"1,0 < T2,l < Ti,i (B7) 

As example of loop configurations that respect this ordering is given in Fig. O As a conse- 
quence, we have: 

To = ri^o and n = r2,i (B8) 

Now we notice that the logarithm of the gauge partition function Zbf,cg{^) in Eq. (1B3P 
contains a sum of differences of the winding angles defined in Eq. (1B6P : 

2Tx\og Zbf,cgW rn ( \ -n ^^,n ^^n ( \ 

TT = [-t>'l,l;2,HT2,lJ - i>'l,l;2,Un,0j + i^l,2;2,2 (T2,1 J " -t>'l,l;2,2 (T2,1 J 

L/\ 

+ ^l,2;2,l(n,o) - -C'l,2;2,l(r2,l) + -Dl,l;2,2(ri,o) " £'l,2;2,2(ri,o)] (B9) 

Further, assuming that the curves Fi and F2 are oriented as in Fig. O if we start from the 
minimum point at Tq = Ti 0, we can isolate in the right hand side of Eq. ( ]B9P the following 
four contributions: 

1. In the time slice ri < i < T2,i the angle which measures the winding of the trajectory 
Fi 1 around the trajectory F2,i is given by the difference -Di,i:2,i(''"2,i) — -Di,i;2,i(ti,o)- 

2. In the region r2,i <t< ti^ only the trajectory Fi continues to evolve, going first up- 
wards with the subtrajectory Fi^i and then downwards with Fi^2- After this evolution, 
the winding angle between the two trajectories Fi and F2 has changed by the quantity 

-Di,2;2,2(''"2,i) — -Dl,l;2,2(''"2,l)- 

3. Next, in the region r2,i > t > Ti^, the winding angle which measures how many times 
the subtrajectory F12 winds up around F2,2 is given by the difference -Di,2;2,i(''"i,o) ~ 

^l,2;2,l(r2,l). 
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4. Finally, in the region ri q > t > T2^o only the second trajectory r2 continues to evolve, 
going first downwards with the curve r2,2 and then upwards with r2,i. The net effect 
of this evolution is that the winding angle between Fi and F2 changes by the quantity 

^l,l;2,2(n,o) - -Di,2;2,2(n,o)- 

It is thus clear that the right hand side of Eq. (1B9I) . apart from a proportionality factor 
iX, counts how many times the trajectory Fi winds around the trajectory F2. If we wish 
to identify the quantity in the right hand side of Eq. flB9p with the Gauss linking number 
x(Fi,F2), we should check for consistency that it takes only integer values as the Gauss 
linking number does. Indeed, it is easy to see that, modulo 27r, the following identities are 
holding: 

^l,l;2,l('''2,l) = -Di^i.2,2(t"2,i) 
-Dl,l;2,2('''l,o) = -Dl,2;2,2('''l,o) 
Di^2;2,2{T2,i) = -Dl,2;2,l(''"2,l) 
i^l,l;2,l(ri,o) = /^l,2;2,l(ri,o) (BIO) 

For example, the first of the above equalities states that the angle formed by the vector 

ri,i — r2,i connecting the subtrajectories Fi 1 and F2,i at the height r2,i is equal to the angle 

formed by the vector ri^i — r2,2 connecting the subtrajectories Fi^i and F2,2 at the same 

height. The reason of this identity is trivial: At that height, the subtrajectories F2,i and 

F2,2 are connected together at the same point. Applying the above relations to Eq. (lB9p . 

one may prove that: 

— ^ = mod 2% (-dH) 

iX 

As a consequence, we can write: 

Zbf,cgW = e^'^(^^'^^) (B12) 

where x(Fi,F2) is the Gauss linking number. Concluding, the above analysis shows that 
also in the Coulomb gauge the BF fields in the polymer partition function (I43p fix the 
topological constraints f ITT]) correctly, in full consistency with the results obtained in the 
covariant gauges. Of course this consistency was expected due to gauge invariance. Yet, it 
is interesting that, using the Coulomb gauge, one may express the Gauss linking number 
invariant in a way that is quite different from the usual form given in Eq. (fTOl) . 
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Appendix C: Prom polymers to anyon field theories 



In this Appendix the passage from the polymer partition function (H3|) to the field the- 
oretical formulation of Eq. ( 146|) is performed. To this purpose, we have to integrate over 
all polymer trajectories r^ /^(ta^/^), a = 1,2 and /q = 1, . . . , 2sa. The standard procedure to 
pass to field theory in polymer physics consists in introducing auxiliary fields. This proce- 
dure works of course also in the present case, but it is complicated by the splitting of the 
trajectories Ta into 2sa subtrajectories. First of all, we have to introduce external sources 
for each subtrajectory as follows: 

J,,,,(x,t)= r'^ dt,jJi^-r,jSta,lJ)6it-ta,lJi-iy^-' (CI) 

Here the coordinates (x, t) are allowed to span the whole M.^ space. Now it is possible to 
write the following identity: 



exp 



- fd^xcfydtdt'JajS^, t)V{x, y)6it - t')Jbj,{y, t'\ 



(C2) 



where a,b = 1,2, la = l,...,2sa, Jb = l,---,2sb. Clearly, the right hand side of the 
above equation can be interpreted as the generating functional of a free scalar field theory 
with propagator ^(x, y; t, t') = V(x, y)(5(t — t'). At this point we notice that the weight 
q-Sev ^]-^at takes into account the short-term interactions in the partition function (H3|) is 
a product of exponents of the kind given in Eq. (1C2I) . Thus, e""^^^ coincides formally with 
the generating functional of a multi-component scalar field theory. The minimum number of 
scalar fields that is necessary to express e""^^^ as a generating functional is 2sa+2sb+2. Let's 
call these fields (/?i^7^(x, t), v92,Jt,(x, t), la = l,...,2sa, Jb = l,...,2sb and 0i(x, t), 02(x, t). 
(pi{'x,t) and 02(x, t) will be responsible for the interaction between monomers belonging to 
different loops, while the 99i/^(x, t)'s and (/J2,Jb(x, t)'s will take into account the interactions 
of monomers belonging to the same loop. 

Remembering that in the present case y(x, y) = Vo(^(x — y), it is possible to verify the 
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validity of the following identity: 



X n / '^^hii^'t) n / ^v'2,j(x,t) 

\J2,Jo = l 



X exp 



X exp 



X exp 



X exp 



2V, 

2s 



2s2 „ 
lJ2 / rf^xdtJ2,j(x,t))(0i(x,t)+V?2,j(x,t)) 

7 T ^ 



(C3) 



where di^j'^, la, I'a = ^, ■ ■ ■ , "^Sa, is the off-diagonal matrix 



o^iai' 



if 4 = I'a 



(C4) 



and (yj\i represents its inverse. Using equation (ICSp . the partition function (H3|) becomes: 



Z(A) = fvB^VC^e-''^^ J 

2si 2s2 

X JJp(^i,,(x,t) JJp(^2,7(x,t) 



V(t)iV(j)2e~^o ^ ^''""^'^'^""'''^'''^^""'''^ 



1=1 



J=l 



X exp 



^ „ / 2si 2s2 

-^ / ^^xdt ^ (/5i,7,(x,t)(^i,7/(x,s)a-^^, + ^ (^2,j2(x,t)^2,4(x,s)a7^^_^, 



2si 1 ^rij(ri,j) 

n / ^ri,,(tM) 

7=1 "'rij(ri,j_i) 
■2*2-1 /.r2..7(r2,.j) 

n / ^r2,j(t2,j; 

J=l >^I"2,/(f"2,J-l) 



ri,2si{-ri,o) 

ri,2si{Tl,2si-l) 
r2,2s2{''"2,o) 

I"2,2s2{''"2,2s2-l) 



^I'l,2si(^l,2si, 



^r2,2s2(^2,252)e" 



'e// 



(C5) 
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with 

2si „^j/ 
1=1 Jti,I-1 



-ly-'gij 



dri 



dti 



+ {-iy-H{(j)2{vi^i{tij)Mj) + V^i,/(ri,/(ti,/),ti,7)) 



t/ri J 






J=l >^^2,J-1 



^y-'92,j 



dT2,. 



dt 



2, J 



+ (-l)-'-'i(0i(r2,j(t2,j),t2,j) + V52,j(r2,j(M.^2,j)) 



+ mC3(r2,j(t2,j),t2,j) + ^:^ ■ C(r2,j(t2,j),t2,j) 



27rdt 



2, J 



(C6) 



Of course, the ends of the trajectories Fi / and T2,.j appearing in the hmits of path integration 
over ri^/(ti_/) and r2,j(t2,j) in Eq. (ICSP are not all independent, because they are subjected 
to the constraints ([3]) and (jlj). We will get rid of these constraints later when passing to the 
field theoretical representation. 

Let us use at this point the so-called complex replica fields defined in Eq. (l47j) . Then, 
the path integrals over the trajectories r^/ may be rewritten as follows 50|: 



2si 1 prij{Tij) 

n / ^ri,,(ti,/) 

7 = 1 >^I"ij(ti,j_i) 

■2S2-1 pr2.j{T2,j) 

n / ^r2,j(t2,j 

J— I >^r2,/(f"2,.j-i) 
■2si-l 



ri,2s]^ (n,o) 

ri,2s]^{Tl,2si-l) 
r2,2s2(''"2,o) 

r2,2s2{''"2,2s2-l) 



^ri,2si(^l,2si) 
:Pr2,2s2(^2,2s2)e"^^^^ 



lim 



1=1 "^ 



X V'i%i(ri,2.i-i(ri,2.i-i),ri,2.i-i)V't2si(ri,i(n,o),T-i,o) 

9 1 



X 



1=1 [^ Jti,I-1 J 

+ i\{-iy-^B^{^, t) + z(02(x, t) + (/Ji,,(x, t)) 






:V-zA(-ir-^B(x,t)) 



2S2-1 



X lim fr /p^2,j:P^;,j^2j(r2,j-i(r2„/-i),r2,j_i)V'2,j(r2,j(r2,j),r2,j) 
J=l -^ 



X ^2*2S2 (r2,2s2-l (^2,2s2-l) , ^2,2s2-l) ^^2,2.2 ("^2,1 (r2,o) , ^2,o) 
2S2 f /•T2,J+1 

X I I exp 



J=2 



IK 



T2,J 



dt{-iy-' /d^xvi/* (x,t) 






+ — (-1)'-'C3(X, t) + 2(01 (X, t) + (^2,/(x, t)) 



^2,j(x,t) 



(C7) 
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Let us note that with the above choice of arguments of the fields if}l*j,ifjij and il'2*j^'^2,j, 
I = 1, . . . , 2si, J = 1, . . . , 2s2, the constraints ([3]) and (jl]) are already taken into account. 
Inserting this result in Eq. ( 1C5I) and integrating out the auxiliary fields ipu, f2,j, 4>i and (p2, 
we obtain the final expression of the polymer partition function given by Eqs. (H6ll . ( HHl) and 

ED. 
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